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ABSTRACT 

We compare the amplitudes of fluctuations probed by the 2dF Galaxy Redshift Survey 
and by the latest measurements of the Cosmic Microwave Background anisotropies. 
By combining the 2dFGRS and CMB data we And the linear-theory rms mass fluc¬ 
tuations in 8 h~^ Mpc spheres to be crsm = 0.73 ± 0.05 (after marginalization over the 
matter density parameter Om and three other free parameters). This normalization is 
lower than the COBE normalization and previous estimates from cluster abundance, 
but it is in agreement with some revised cluster abundance determinations. We 
also estimate the scale-independent bias parameter of present-epoch Lg = 1.9L* 
APM-selected galaxies to be b{Ls,z = 0) = 1.10 ± 0.08 on comoving scales of 
0.02 < k < 0.15/iMpc“^. If luminosity segregation operates on these scales, 
galaxies would be almost un-biased, b{L^, z = Q) ~ 0.96. These results are derived by 
assuming a flat ACDM Universe, and by marginalizing over other free parameters and 
fixing the spectral index n = \ and the optical depth due to reionization r = 0. We 
also study the best fit pair (Um, &), and the robustness of the results to varying n and 
T. Various modelling corrections can each change the resulting b by 5-15 per cent. 
The results are compared with other independent measurements from the 2dFGRS 
itself, and from the SDSS, cluster abundance and cosmic shear. 
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1 INTRODUCTION 


The 2dF Galaxy Redshift Survey (2dFGRS) has now mea¬ 
sured over 210 000 galaxy redshifts and is the largest existing 
galaxy redshift survey (Colless et al. 2001). A sample of this 
size allows large-scale structure statistics to be measured 
with very small random errors. Two other 2dFGRS papers, 
Percival et al. (2001; hereafter POl) and Efstathiou et al. 
(2001; hereafter E02) have mainly compared the shape of the 
2dFGRS and GMB power spectra, and concluded that they 
are consistent with each other (see also Tegmark, Hamil¬ 
ton & Xu 2001). Here we estimate the amplitudes of the rms 
fluctuations in mass usm and in galaxies crsg. More precisely, 
we consider the ratio of galaxy to matter power spectra, and 
use the ratio of these to define the bias parameter: 


_ Pggj^) 

Pruru (k) ’ 


( 1 ) 


As defined here, b is in principle a function of scale. In prac¬ 
tice, we will measure the average value over the range of 
wavenumbers 0.02 < k < 0.15/iMpc“^. On these scales, 
the fluctuations are close to the linear regime, and there 
are good reasons (e.g. Benson et al. 2000) to expect that b 
should tend to a constant. In this study, we will not test the 
assumption that the biasing is scale-independent, but we do 
allow it to be function of luminosity and redshift. 

A simultaneous analysis of the constraints placed on 
cosmological parameters by different kinds of data is essen¬ 
tial because each probe - e.g. GMB, Type la supernovae 
(SNe la), redshift surveys, cluster abundance, and pecu¬ 
liar velocities - typically constrains a different combination 
of parameters (e.g. Bahcall et al. 1999; Bridle et al. 2000, 
2001a; E02). A particular case of joint analysis is that of 
galaxy redshift surveys and the Gosmic Microwave Back¬ 
ground (GMB). While the GMB probes the fluctuations in 
matter, the galaxy redshift surveys measure the perturba¬ 
tions in the light distribution of particular tracer (e.g. galax¬ 
ies of certain type). Therefore, for a fixed set of cosmological 
parameters, a combination of the two can tell us about the 
way galaxies are ‘biased’ relative to the mass fluctuations 
(e.g. Webster et al. 1998). 

A well-known problem in estimating cosmological pa¬ 
rameters is the degeneracy of parameters, and the choice of 
free parameters. Here we consider three classes of parame¬ 
ters: 

(i) Parameters that are fixed by theoretical assumptions 
or prejudice (which may be supported by observational evi¬ 
dence). Here we assume a flat Universe (i.e. zero curvature), 
and no tensor component in the GMB (for discussion of the 
degeneracy with respect to these parameters see E02). 

(ii) The ‘free parameters’ that are of interest to address 

a particular question. For the joint 2dFGRS & GMB analy¬ 
sis presented here we consider five free parameters: the mat¬ 
ter density parameter Hm, the linear-theory amplitude of 
the mass fluctuations usm, the present-epoch linear biasing 
parameter b{Ls ,z = 0) (for the survey effective luminosity 
La ~ 1.9Z/*), the Hubble constant h = Ho/(100kms“^), 
and the baryon density parameter ujh = As we are 

mainly interested in combinations of asm, b and Hm, we 
shall marginalize over the remaining parameters. 

(iii) The robustness of the results to some ‘extra pa¬ 
rameters’, that are uncertain. Here we consider the optical 


depth r due to reionization (see below) and the primordial 
spectral index n. We use as our canonical values r = 0 and 
n = 1, but we also quote the results for other possibly real¬ 
istic values, r = (0.05, 0.2) and n = (0.9, 1.1). 

The outline of this paper is as follows. In Section 2 
we derive ngg from the 2dFGRS alone, taking into ac¬ 
count corrections for redshift-space distortion and for epoch- 
dependent and luminosity-dependent biasing. In Section 3 
we derive agm from the latest GMB data. In Section 4 we 
present a joint analysis of 2dFGRS & GMB. Finally, in Sec¬ 
tion 5 we compare and contrast our measurements with 
other cosmic probes. 


2 THE AMPLITUDE OF THE 2dFGRS 

FLUCTUATIONS 

2.1 crfg from the fitted power spectrum 

An initial estimate of the convolved, redshift-space power 
spectrum of the 2dFGRS has already been determined 
(POl), using the Fourier-transform-based technique de¬ 
scribed by Feldman, Kaiser & Peacock (1994; hereafter 
FKP) for a sample of 160 000 redshifts. On scales 0.02 < k < 
0.15 h Mpc”*^, the data are robust and the shape of the power 
spectrum is not affected by redshift-space or non-linear ef¬ 
fects, though the amplitude is increased by redshift-space 
distortions (see later). We use the resulting power spectrum 
from POl in this paper to constrain the amplitude of the 
fluctuations. 

As explained above, we define the bias parameter as the 
square root of the ratio of the galaxy and mass power spectra 
on large scales. We shall assume that the mass power spec¬ 
trum can be described by a member of the family of models 
dominated by Gold Dark Matter (GDM). Such models tra¬ 
ditionally have their normalization described by the linear- 
theory value of the rms fractional fluctuations in density 
averaged in spheres of 8/i“^Mpc radius: crgm. It is there¬ 
fore convenient to define a corresponding measure for the 
galaxies, ngg, such that we can express the bias parameter 
as 

b=^. ( 2 ) 

ngm 

The scale of 8 Mpc was chosen historically because ngg ~ 
1 from the optically-selected Lick counts (Peebles 1980), so it 
may seem impossible by definition to produce a linear-theory 
erg for galaxies. In practice, we define ergg to be the value 
required to fit a CDM model to the power-spectrum data 
on linear scales (0.02 < k < 0.15 h,Mpc“^). From this point 
of view, one might equally well specify the normalization via 
e.g. CT 2 o; however, the ag parameter is more familiar in the 
context of CDM models. The regions of the power spectrum 
that generate the erg signal are at only slightly higher k than 
our maximum value, so no significant uncertainty arises from 
extrapolation. A final necessary complication of the notation 
is that we need to distinguish between the apparent values of 
agg as measured in redshift space (erfg) and the real-space 
value that would be measured in the absence of redshift- 
space distortions (cr|^). It is the latter value that is required 
in order to estimate the bias. 

The 2dFGRS power spectrum (Fig. 1) is fitted in POl 
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Figure 1. The observed (i.e. convolved with the window function) 
2dFGRS power spectrum (as derived in POl). The solid line shows 
a linear theory ACDM fit (convolved with the window function) 
with Q,uih = 0.2,Oi,/^m = 0.15, h = 0.7, n = 1 and best-fitting 
a^^{zs, Ls) = 0.94. Only the range 0.02 < k < 0.15hMpc“^ 
is used at the present analysis (roughly corresponding to CMB 
harmonics 200 < i < 1500 in a flat Qm = 0.3 Universe). The good 
fit of the linear theory power spectrum at fc > 0.15 /iMpc“^ is due 
to a conspiracy between the non-linear growth and finger-of-god 
smearing (integrating over the observed P{k) therefore provides 
another way of estimating the normalization, giving ~ 0.95). 



Figure 2. The likelihood function of 2dFGRS as a function of 
of the galaxy fluctuation amplitude in redshift space {La, Za) 
and the present epoch rim. The marginalization over the Hubble 
constant is done with a Gaussian centred at h = 0.7 and standard 
deviation of 0.07. Other parameters are held fixed (n = l,ti;b = 
0.02). The contours contain 68 per cent and 95 per cent of the 
probability. 

over the above range in k, assuming scale-invariant primor¬ 
dial fluctuations and a ACDM cosmology, for four free pa¬ 
rameters: 0,mh, Ob/Dm, h and the redshift space afg (us¬ 
ing the transfer function fitting formulae of Eisenstein & 
Hu 1998) . Assuming a Gaussian prior on the Hubble con¬ 
stant h = 0.7 db 0.07 (based on Freedman et al. 2001) the 
shape of the recovered spectrum within the above fc-range 
was used to yield 68 per cent confidence limits on the shape 


parameter Qmh = 0.20 ± 0.03, and the baryon fraction 
Db/S2m = 0.15 ± 0.07, in accordance with the popular ‘con¬ 
cordance’ modelj^. Although the ACDM model with compa¬ 
rable amounts of dark matter and dark energy is rather eso¬ 
teric, it is remarkable that the 2dFGRS measurement shows 
such good consistency with other cosmological probes such 
as CMB, SNe, and Big Bang Nucleosynthesis (BBN). 

We And that afg depends only weakly on the other 
three parameters, with the strongest correlation with flm- 
For fixed ‘concordance model’ parameters n = l,Dm = 
1 — Qa ~ 0.3, uih = 0.02 and a Hubble constant h = 0.70, 
we find that the amplitude of 2dFGRS galaxies in redshift 
space is G^^{La,Za) = 0.94 (when all other parameters are 
held fixed, the formal errors are unrealistically tiny, only a 
few percent, and hence we do not quote them). In the FKP 
method, the normalization of the power spectrum depends 
on the radial number density and weighting function, and 
a number of different methods have been suggested for cal¬ 
culating the normalization (Sutherland et al. 1999) using 
a random catalogue designed to Poisson sample the survey 
region. POl tried all of the suggested methods for the 2dF- 
GRS data and found no significant change in the power spec¬ 
trum normalization. Therefore, although this calculation re¬ 
mains a potential cause of systematic error in the power- 
spectrum normalization, we shall assume hereafter that the 
main uncertainty in the derived bias derives from the uncer¬ 
tain cosmological model that is needed in order to connect 
the galaxy power spectrum with the mass power spectrum 
from the CMB. 

On keeping flm free and marginalizing over h with a 
Gaussian prior (with h = 0.7± 10 per cent) we obtain Fig. ^ 
The external constraint on h that we impose translates to 
a constraint on flm through the 2dFGRS sensitivity to the 
matter power spectrum shape, which is roughly flmh. On 
marginalizing over flm we find crfg = 0.94 ± 0.02, in agree¬ 
ment with the best-fit non-marginalized resultQ 

2.2 Corrections for redshift and luminosity effects 

In reality, the effective redshift for the POl analysis is not 
zero, but Za ~ 0.17. This is higher than the median redshift 
of 2dFGRS {zm ~ 0.11) due to the weighting scheme used 
in estimating the power spectrum. Similarly, La — 1.91/*, 
rather than the La — A* that would apply for a flux-limited 
sample. We can then derive ag^{La, Za) directly, but for com¬ 
parison with other studies we make further steps of calcu- 

^ As shown in POl, the likelihood analysis gives a second (non¬ 
standard) solution, with Qmh ~ 0.6, and the baryon fraction 
flb/flm = 0.4, which generates baryonic ‘wiggles’. We ignore this 
case at the present analysis and using the likelihood function over 
the range 0.1 < flmh < 0.3, 0.0 < Ob/flm < 0.4, 0.4 < h < 0.9 
and 0.75 < Ugg < 1.14. We also note that even if there are features 
in the primordial power spectrum, they would get washed out by 
the 2dFGRS window function (Elgaroy, Gramann & Lahav 2002) 
^ We emphasise again that here agg is the linear-theory nor¬ 
malization, not the observed non-linear cggNL- For example, the 
2dFGRS correlation function of Norberg et al. (2001a) can be 
translated to a non-linear = 0.87±0.07, at an effective 

redshift of approximately 0.07. In practice, nonlinear corrections 
to (Tg are expected to be relatively small for CDM-like spectra 
(see Fig. 1). 
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Figure 3. A compilation of the latest CMB data points from 
COBE, Boomerang, Maxima and DASI against spherical har¬ 
monic £. The line shows the predicted angular power spectrum 
for a ACDM model with n = 1, = 1 — Ha = 0.3, cJb = 0.02 

(BBN value), h = 0.70, r = 0.0, and the best-fit normalization to 
the given CMB data points (Tgm = 0.83. Note that this normaliza¬ 
tion is lower than the traditional COBE-only normalization (see 
Table 1). A similar model is also the best fit to the shape of the 
2dFGRS galaxy power spectrum (Fig. 1). 

lating cr|^(Ls, z = 0) and then cr||(L*, z = 0). This requires 
corrections which depend on the nature of galaxy forma¬ 
tion and on clustering with redshift. Some of the corrections 
themselves depend on cosmological parameters, and our pro¬ 
cedure solves for the best-fitting values in a self-consistent 
way. 

We start by evaluating the conversion from redshift 
space to real space at the survey effective redshift for 
galaxies with effective luminosity Ls'. 

ag{Ls,Zs) = ai^{Ls,Zs)/K^^^[l3{L,,z,)] . (3) 

where 

7f[/3] = l + |/3+i/3' (4) 

is Kaiser’s (1987) factor, derived in linear theory and the 
distant-observer approximatiorf]. The dependence of f3 on 
redshift can be written as: 

P{Ls,Zs) ~nJ’-\zs)/b{Ls,Zs), (5) 

assuming linear biasing [for more general biasing schemes 
see e.g. Dekel & Lahav (1999) and references therein]. 

The evolution of the matter density parameter with red¬ 
shift is 

Qm{z) = nu,{l + zf {H/Ho)-^ (6) 

with 

^ More precisely, the redshift-space distortion factor depends on 
the auto power spectra Pmm{h) and Pgg(fc) for the mass and the 
galaxies, and on the mass—galaxies cross power spectrum Pmg{k) 
(Dekel Sz Lahav 1999; Pen 1998; Tegmark et al. 2001). The model 
of equations (3—5) is only valid for a scale-independent bias factor 
b that obeys Pgg(fc) = bPmg(h) = 5^Pmm(fc)- 


{H/Hof = [nm(l + 2 )® -H (1 - - flA)(l + zf + flA] .(7) 

The variation of b{z) = (Jsg{z)/asiniz) with redshift is 
even more model-dependent. We assume that the mass fluc¬ 
tuations grow as (T 8 m( 2 ) = cr8m(0)D(2), where D{z) (nor¬ 
malized to 1 at 2 = 0) is the growing mode of fluctuations 
in linear theory [it depends on flm and Ha, e.g. Peebles 
(1980)]. 

We also assume that galaxy clustering weakly evolves 
over 0 < z < 0.2, i.e. asg{Ls,0) ~ asg{Ls,z). We shall 
refer to this simple model as the ‘constant galaxy clustering 
(CGC) model’. Simulations suggest (e.g. Kauffmann et al. 
1999; Blanton et al. 2000; Benson et al. 2000; Somerville 
et al. 2001) that even if the clustering of dark matter halos 
evolves slightly over this range of redshifts, galaxy clustering 
evolves much less. Indeed, observationally there is only a 
weak evolution of clustering of the overall galaxy population 
over the redshift range 0.1 < 2 < 0.5 (e.g. the CNOC2 
survey: Shepherd et al. 2001). Therefore in our simple CGC 
model (for any luminosity): 

biLs,Zs) = b{Ls,0)/Dizs). (8) 

There are of course other possible models for the evolu¬ 
tion of galaxy clustering with redshift, e.g. the galaxy con¬ 
serving model (Fry 1996). This model describes the evolu¬ 
tion of bias for test particles by assuming that they follow 
the cosmic flow. It can be written as: 

b{Ls, zs) = l + [fe(Ls, 0) - 1]/D{zs). (9) 

More elaborate models exist, such as those based on a merg¬ 
ing model (e.g. Mo & White 1996; Matarrese et al. 1997; 
Magliocchetti et al. 1999) or numerical and semi-analytic 
models (Benson et al. 2000; Somerville et al. 2001). 

To estimate the magnitude of these effects we consider 
the 2dFGRS effective redshift Za = 0.17. For a Universe with 
present-epoch flm = 0.3 and Ha = 0.7 we get nin(2s) = 0.41 
and D{za) = 0.916, hence for the CGC model with = 
0.3, b{La,Za) = 1.09 b{La,0) and P{La, Za) = 1.10 P{La, 0). 

On the other hand, we can also relate the amplitude 
of galaxy clustering to the present-epoch mass fluctuations 
(Tsm, which can be estimated from the CMB (see below): 

agg{La,0) = b{La,0)asni{0) . (10) 

Hence by combining equations (3), (8) and (10) we can solve 
for b{La, 0). 

Finally, there is the issue of luminosity-dependent bi¬ 
asing. Although controversial for some while, this effect has 
now been precisely measured by the 2dFGRS (Norberg et al. 
2001a,2002); see also recent results from the SDSS (Zehavi 
et al. 2001). Norberg et al. (2001a) found from correlation- 
function analysis that on scales ^ 10 h~^ Mpc 

6(L,0)/6(L.,0) = 0.85-f 0.15(L/L*) . (11) 

If we assume that this relation also applies in the lin¬ 
ear regime probed by our P{k) on scales 0.02 < k < 
0.15/iMpc”*^, then the linear biasing factor for L* galax¬ 
ies at redshift zero is 1.14 smaller then that for the 2dFGRS 
galaxies with effective survey luminosity La = 1.9L* {K- 
corrected). However, this is a source of uncertainty, and ul¬ 
timately it can be answered with the complete 2dFGRS and 
SDSS surveys by calculating the power spectra in luminosity 
bins. 
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The luminosities in equation (11) have been K- 
corrected and also corrected for passive evolution of the 
stellar populations, but the clustering has been measured at 
various median redshifts (for galaxies at the redshift range 
0.02 < 2 : < 0.28). Possible variation of galaxy clustering with 
redshift is still within measurement errors of Norberg et al. 
(2001a,2002). For simplicity we shall assume in accord with 
our CGC model that this relation is redshift-independent 
over the redshift range of 2dFGRS (z<0.2). We see that 
the effects of redshift-space distortion and luminosity bias 
are quite signihcant, at the level of more than 10 per cent 
each. 


3 THE CMB DATA 


The GMB fluctuations are commonly represented by the 
spherical harmonics Ci. The connection between the har¬ 
monic £ and k is roughly 


I ~ kdA 


( 12 ) 


where for a flat Universe the angular distance to the last 
scattering surface is well approximated by (Vittorio & Silk 
1991): 


dA 


2 c 

Ho nr ■ 


(13) 


For rim = 0.3 the 2dFGRS range 0.02 < k < 0.15 hMpc”^ 
corresponds approximately to 200 < £ < 1500, which is 
well covered by the recent GMB experiments. We obtain 
theoretical GMB power spectra using the GMBFAST code 
(Seljak & Zaldarriaga 1996). 

The latest GMB measurements from Boomerang (Net- 
terfield et al. 2001, de Bernardis et al. 2002), Maxima (Lee 
et al. 2001; Stomper et al. 2001) and DASI (Halverson et 
al. 2002; Pryke et al. 2002) suggest three acoustic peaks. 
Parameter htting to a AGDM model indicates consistency 
between the different experiments, and a best-ht Universe 
with zero curvature, and an initial spectrum with spectral 
index n ~ 1 (e.g. Wang, Tegmark & Zaldriaga 2001, here¬ 
after WTZ; E02 and references therein). Unlike the earlier 
Boomerang and Maxima results, the new data also show 
that the baryon contribution is consistent with the Big Bang 
Nucleosynthesis value uih — 0.02 (O’Meara et al. 2001). 

Various GMB data sets can be combined in different 
ways (e.g. Jaffe et al. 2001; Lahav et al. 2000). Here we 
consider two compilations of GMB data: 

(a) a compilation of GOBE (8 points). Boomerang, 
Maxima and DASI (hereafter GBMD). The total number 
of data points in this compilation is 49, plotted in Fig. 

(b) a compilation of 24 AT/T data points from Wang et 
al. (2001; WTZ), which is based on 105 band-power measure¬ 
ments of almost all available GMB experiments (including 
the latest Boomerang, Maxima, and DASI data). 

Both compilations take into account the calibration er¬ 
rors, which are crucial for estimating the amplitude of fluc¬ 
tuations. For our compilation (a) we use a fast method 
for marginalization over calibration and beam uncertain¬ 
ties that assumes a Gaussian prior on the calibration and 
beam corrections (Bridle et al. 2001b). We apply the usual 


Table 1. Normalizations for matter fluctuations derived by fitting 
to CMB data alone. In all the entries (unless otherwise stated) 
other parameters are fixed at Um = 1 — Ga = 0.3, wt, = 0.02, n = 
1 and h = 0.7. The first three entries were derived via CMBFAST 
using the COBE points according to the normalization procedure 
of Bunn &; White (1996). The other entries were derived by the 
best fit multi-variate (including the covariance matrix and 
window functions) for the WTZ and CBMD data points (the 
goodness of fit is e.g. x^ = 35 for 24 points in the fourth entry). 
Quoted error bars are 1-sigma. In some cases formal errors are 
not quoted as they are unrealistically tiny (few percent) when 
other parameters are held fixed. Note that both the WTZ and 
CBMD compilations give normalization lower than the COBE- 
only normalization. 


Data 


COBE (r = 0), 

0.90 

COBE (r = 0.05), 

0.93 

COBE (r = 0.20), 

0.98 

WTZ (r = 0), 

0.77 

WTZ (r = 0.05), 

0.80 

WTZ (r = 0.20), 

0.92 

CBMD (r = 0) 

0.83 

CBMD(t = 0, marg. over h = 0.7 ± 0.07 ) 

0.71 ± 0.07 

CBMD (r = 0, marg. over h = 0.7 zb 0.07 & 11m) 

0.68 ± 0.07 

CBMD (r = 0) + 2dF, marg. over h, Sim, uJb ^ b 

0.73 ± 0.05 


multi-variate procedure (e.g. Hancock et al. 1998), tak¬ 
ing into account the window functions and the covariance 
matrix (when available). Since the Boomerang and Maxima 
window functions and correlation matrices are not yet avail¬ 
able, we assume that the data points are uncorrelated and 
use top-hat window functions (as did WTZ). This assump¬ 
tion is validated by the fact that we obtain sensible values 
of for the best-fitting models. 


3.1 CMB-only fits 

We first consider the constraints arising from the GMB data 
alone. Table 1 summarizes various estimates for crsm from 
the above two new data sets. Note that these differ from the 
normalization returned by GMBFAST when only the GOBE 
data are considered. This Table also illustrates the sensitiv¬ 
ity of the results to the optical depth to reionization r (see 
below). We see that differences in data sets and in assump¬ 
tions on other parameters can easily lead to uncertainties of 
~ 10 per cent in the resulting asm- Note that the normaliza¬ 
tions derived from WTZ are lower than those derived from 
our compilation. This reflects that fact that WTZ chose to 
adjust downwards the calibrations of the principal datasets 
that we prefer to adopt. We incorporate the calibration un¬ 
certainties, but make no such adjustment. 

Fig. ^ (dashed lines) shows the likelihood as a function 
of (Qm, nsm) after marginalization over the Hubble constant 
is done with a Gaussian with h = 0.7 ± 0.07, while keep¬ 
ing other parameters fixed (n = l,a;b = 0.02, r = 0.0). We 
note that for a fixed Qm = 0.3 on this diagram the result¬ 
ing CTsm ~ 0.7 is lower that the value we obtained above 
(crsm ~ 0.8) when fixing the Hubble constant h = 0.7 and 
other parameters. This illustrates the sensitivity of the re- 
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Figure 4. The likelihood function of CMB alone (dashed lines) 
in terms of the mass fluctuation amplitude agm and the present 
epoch rim- The marginalization over the Hubble constant is done 
with a Gaussian centred at h = 0.7 and standard deviation of 
0.07. Other parameters are held fixed (n = 1, lUn = 0.02, r = 0.0). 
The contours are for (two-parameter) 68 per cent and 95 per cent 
confidence intervals. The solid lines show the contours (68 per 
cent and 95 per cent) for the joint 2dFGRS&CMB analysis, after 
marginalization over h, b(Ls,0) and lUu- Other parameters are 
held fixed (n = 1,t = 0.0). Note that the contours of 2dFGRS & 
GMB are much tighter than when using CMB alone. Two recent 
extreme cluster abundance determinations are overlayed as the 
upper dotted line (Pierpaoli et al. 2001) and the lower dotted line 
(Viana et al. 2002). 

suits from the CMB alone to the Hubble constant. The ex¬ 
ternal constraint on h we have imposed cuts off the contours 
at low and high flm. This is due to the constraint on 
that exists from CMB data: a constraint on h thus translates 
to a constraint on flm. Completing the marginalization over 
flm we find crsm = 0.68 ± 0.07. Note that since we assume 
that the Universe is flat, there are additional constraints on 
our free parameters that come from the position of the first 
acoustic peak that make our error bars slightly smaller than 
studies that marginalize over the curvature of the Universe 
as well. We overlay in Fig. ^ the constraints from cluster 
abundance obtained recently by various authors. The clus¬ 
ter abundance constraint is fortunately orthogonal to the 
CMB constraint, but the spread in normalization values is 
quite large. It is interesting that some of the latest estimates 
are in good agreement with our estimates from the CMB and 
2dFGRS-l-CMB (see further discussion below). 

4 COMBINING 2dFGRS & GMB 

When combining 2dFGRS and GMB data the parameteri- 
sation for the log-likelihoods is in five parameters: 

ln£tot = ln£ 2 dFGRs[f^m,h, Wb,Cr8m,6(Ts,0)] 

+ In £cmb [fl m 5 ; *^ 8 m] 7 

where £ 2 dFGRS and £gmb are the likelihood functions for 
2dFGRS and the CMB. 

The 2dFGRS likelihood function takes into account the 
redshift-space distortions, the CGC biasing scheme, and the 
redshift evolution of f2in(«). Here we use our compilation of 


49 CMB data points (shown in Fig. 3). Other parameters 
are held fixed (n = 1 , t = 0 ). 

Fig. I (solid lines) shows the 2dFGRS-bCMB likeli¬ 
hood as a function of (flinjUsm), after marginalization over 
h, b[Ls, 0) and cuh- The peak of the distribution is consistent 
with the result for the CMB alone (shown by the dashed lines 
in Fig. but we see that the contours are tighter due to the 
addition of the 2dFGRS data. Further marginalization over 
Sim gives (Tsm = 0.73 db 0.05. The importance of adding the 
shape information of the 2dFGRS power spectrum is that it 
requires no external prior for h and uJh, unlike deriving usm 
from CMB alone (Table 1). Our result is very similar to the 
value (Tsm — 0.72 derived in E02 using the WTZ data set 
and after marginalizing over the raw 2dFGRS amplitude of 
the power spectrum and other parameters. 

To study the biasing parameter we marginalize the 2dF- 
GRS likelihood over h,uih and agm- Other parameters are 
held fixed (n = 1, r = 0). The resulting likelihood as a func¬ 
tion of [Sim, 5(Ts, 0)] is shown (by solid contours) in Fig. ^ 
Further marginalizing over Sim gives b{Ls,0) — 1.10 ± 0.08 
(1-sigma). With Fry’s biasing scheme (equation 9) b(Ls,0) 
is increased by 8 per cent. 

The effect of changing the spectral index to n = 0.9 is 
shown (by dashed lines) in Fig. H (with r = 0). Results for 
n = 0.9 and n = 1.1 with further marginalization over Sim 
are given in Table 2, showing that b{Ls,0) is slightly down 
and up respectively relative to the standard n = 1 case. 
We see that when we fit CMB data over a wide range of £, 
the effect of changing n is small. This is in contrast with the 
large variation of fitting the normalization with COBE only, 
where for the concordance model usm = (0.72, 0.90,1.13) for 
n = (0.9,1.0,1.1), respectively. 

We also tested sensitivity to the optical depth r. Re¬ 
cent important constraints come from the spectra of SDSS 
quasars, suggesting r > 0.03 — 0.04 (Becker et al. 2001; Fan 
et al. 2002). For fixed n = l,Wb = 0.02, and marginaliza¬ 
tion over SlmjUsm and h we get b{La,0) = 1.06 ± 0.09 for 
T = 0.05, i.e. lower by 4 per cent compared with the case 
of r = 0.0. Note that setting Wb = 0.02 or marginalizing 
over it makes little difference to b{Ls,0). The effect of the 
optical depth is indeed expected to increase agm by a factor 
exp(T), and hence to decrease b by that factor, about 5 per 
cent in the case of r = 0.05 (corresponding to redshift of 
reionization Zr — 8 for the concordance model parameters; 
e.g. Griffiths & Fiddle 2001). 

Other possible extra physical parameters may also 
slightly affect our result. For example, a neutrino with mass 
of 0.1 eV (e.g. Hu, Eisenstein & Tegmark 2001; Gawiser 
2001 ) would reduce agm by a few per cent. 

Finally, to translate the biasing parameter from Ls to 
e.g. L* galaxies one can either assume (somewhat ad-hoc) 
no luminosity segregation on large scales, or divide by the 
factor 1.14 (equation 11) that applies on small scales. E.g. 
using the fully marginalized result b(Ls,0) ~ 1.10 we get 
6(1/*,0) ~ 0.96, i.e. a slight anti-bias. Overall, our results 
can be described by the following formula: 

6 (L*, 2 = 0) = (0.96 ± 0.08) exp[—r -|- 0.5(n — 1)]. (15) 
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Figure 5. The result of a joint likelihood 2dFGRS+CMB (solid 
lines). The marginalization (without any external priors) is over 
h,u)\^ and (Tgm- Other parameters are held fixed (n = 1 ,t = 0). 
The contours are for (two-parameter) 68 per cent and 95 per cent 
confidence intervals. The dotted lines represent the 1-sigma enve¬ 
lope for /9(Ls,0), based on /3(Ls,Zs) = 0.43 ± 0.07 from Peacock 
et al (2000) and the CGC model. The result of a joint likelihood 
2dFGRS-|-CMB for n = 0.9 is marked by the dashed contours (68 
per cent and 95 per cent). 


Table 2. The biasing parameter b{Ls, 2 = 0) from the full Maxi¬ 
mum Likelihood solution (equation 14), and marginalization over 
[h, tab) u-gm, flm) without any external priors (apart from the third 
entry, where tab = 0.02). 


Data 

2dFGRS-l-CBMD (n = 1.0, r = 0) 
2dFGRS-l-CBMD (n = 0.9, r = 0) 
2dFGRS-l-CBMD (n = 1.1, r = 0) 
2dFGRS-l-CBMD (n = 1.0, r = 0.05, tab = 0.02) 


b(Ls,0) 

1.10 ±0.08 
1.08 ±0.09 
1.15 ±0.09 
1.06 ±0.09 


5 COMPARISON WITH OTHER 

MEASUREMENTS 

5.1 Other estimates of 2dFGRS amplitude of 
fluctuations 

An independent measurement from 2dFGRS comes from 
redshift-space distortions on scales < 10fi“^Mpc (Peacock 
et al. 2000). This gives /3(Z/s,2s) = 0.43±0.07. In Fig. |^we 
show this constraint, after translating it to f3{Ls, 2 = 0) via 
the CGC model. We see consistency with our present analy¬ 
sis at the level of 1-sigma. Using the full likelihood function 
in the (6, flm) plane (Fig. ^ ) we derive a slightly larger (but 
consistent) value, /3(Ls,2s) — 0.48 ± 0.06. 

A study of the bispectrum of the 2dFGRS (Verde et 
al. 2001) on smaller scales (0.1 < k < 0.5 h Mpc“^) sets 
constraints on deviations from linear biasing, and it gives a 
best-fit solution consistent with linear biasing of unity. The 
agreement with the result of the present paper is impressive, 
given that the methods used are entirely different. In fact, 
by matching the two results one can get constraints on e.g. 
r <0.2. 


5.2 Comparison with other independent 
measurements 

5.2.1 SDSS 

Maximum Likelihood analysis of the early Sloan Digital 
Sky Survey (SDSS) by Szalay et al. (2001) finds from the 
projected distribution of galaxies in the magnitude bin 
20 < r* < 21 (median redshift Zm = 0.33) a shape pa¬ 
rameter F ~ ilmh = 0.183 ± 0.04 and a linear real-space 
cr |^(2 = 0) = 0.785 ± 0.053 (1-sigma errors), assuming a flat 
ACDM model with flm = 1 — Da = 0.3, for the case of no 
evolution of galaxy clustering, equivalent to our CGC model 
(see also Dodelson et al. 2001). To convert the SDSS r* mag¬ 
nitude we use models similar to those given in Norberg et 
al. (2001b), where we find that for Zm — 0.33, 6j ~ r* ± 1 
(for the mix of galaxy populations). Hence at that redshift 
r* = 20 corresponds to absolute M(,j ~ —19.4 (in a flat 
Dm = 0.3 Universe), which with appropriate K and evo¬ 
lution correction gives a rest-frame Mbj — —19.6. This is 
in fact very close to L* of the 2dFGRS. Hence the derived 
SDSS value, (j|^(L *,2 = 0) = 0.785, is in accord with the 
real-space values of we get from 2dFGRS. 


5.2.2 Cluster abundance 

A popular method for constraining asm and Dm on scales of 
~ 10 Mpc is based on the number density of rich galaxy 
clusters. Four recent analyses span a wide range of values, 
but interestingly they are all orthogonal to our CMB and 
2dF constraints (Fig. 4) . 

Pierpaoli, Scott & White (2001) derived a high value, 
while Seljak (2001), Reiprich & Boehringer (2002), and 
Viana, Nichol & Liddle (2002) found lower values: 


asm 0.50 Dm"° ® 

asm 0.44 Dm"°'^'‘ 

asm 0.43 Dm"°'®® 

n oo o —0.48+0.27r2m 
<78m — U.oO i im 


(16) 


respectively. For Dm = 0.3 these results correspond to 
(Jsm — 1.02; 0.75; 0.68; 0.61 respectively (with typical er¬ 
rors of 10 per cent). The high value agrees with numerous 
earlier studies by Eke et al. (1998) and others which were 
based on temperature functions, and it remains to be under¬ 
stood why the recent values are so low. The discrepancy be¬ 
tween the different estimates is in part due to differences in 
the assumed mass-temperature relation. The cluster physics 
still needs to be better understood before we can conclude 
which of the above results is more plausible. We see in Fig. 
4 that the lower cluster abundance results are actually in 
good agreement with our value from the 2dFGRS±CMB, 
asm ^ 0.73 ± 0.05. 


5.2.3 Cosmic shear 

The measurements of weak gravitational leasing (cosmic 
shear) are sensitive to the amplitude of the matter power 
spectrum on mildly non-linear scales. Van Waerbeke et al. 
(2001), Rhodes, Refregier & Groth (2001) and and Bacon et 
al. (2002) find respectively 
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(Tg^ -0.43 ® 

(Tg^ (17) 

(Jgn, -0.43 

(with errors of about 20 per cent). These estimates are 
higher than the ngm value that we obtain from 2dF- 
GRS+CMB, but note the large error bars in this recently 
developed method. 


6 DISCUSSION 

We have combined in this paper the latest 2dFGRS and 
GMB data. The first main result of this joint analysis is the 
normalization of the mass fluctuations, ngm = 0.73 ± 0.05. 
This normalization is lower than the COBE normalization 
and previous estimates from cluster abundance, but it is ac¬ 
tually in agreement with recently revised cluster abundance 
normalization. The results from cosmic shear are still some¬ 
what higher, but with larger error bars. 

The second result is for the biasing parameter for 
optically-selected Ls galaxies, b(Ls,0) = 1.10 ±0.08, which 
is just consistent with no biasing (‘light traces mass’) on 
scales of tens of Mpc. When translated to L* via a cor¬ 
rection valid for small scales we get a slight anti-bias, 
5(1/*,0) — 0.96. Although biasing was commonly neglected 
until the early 1980s, it has become evident that on scales 
<10/i“^Mpc different galaxy populations exhibit different 
clustering amplitudes, the so-called morphology-density re¬ 
lation (e.g. Dressier 1980; Hermit et al. 1996; Norberg et al. 
2002). Biasing on small scales is also predicted in the simula¬ 
tions of hierarchical clustering from CDM initial conditions 
(e.g. Benson et al. 2000). It is important therefore to pay at¬ 
tention to the scale on which biasing operates. Our result of 
linear biasing of unity on scales ( > 10 Mpc) is actually 
in agreement with predictions of simulations (e.g. Blanton 
et al. 2000 Benson et al. 2000; Somerville et al. 2001). It 
was also demonstrated by Fry (1996) that even if biasing 
was larger than unity at high redshift, it would converge 
towards unity at late epochs (see equation 9). 

We note that in deriving these results from the 2dFGRS 
and the GMB, we have had to consider various corrections 
due to astrophysical and cosmological effects: 

• redshift-space distortions cause the amplitude in red- 
shift space to be ~ 15 per cent larger than that in real 
space. 

• The evolution of biasing with redshift (for our simple 
constant galaxy clustering model) gives a biasing that is 
~ 10 per cent higher at Zs = 0.17 than at redshift zero. 

• If luminosity-dependent biasing also holds on large 
scales then the biasing parameter b(Ls = 1.9Z/*) is — 15 
per cent higher than that of L* galaxies. 

• On the GMB side, an optical depth r = 0.05 due to 
reionization reduces the derived biasing parameter fe by — 5 
per cent. Ghanging the spectral index from n = 1 to n = 0.9 
(for both the GMB and 2dFGRS) also reduces 6 by ~ 5 per 
cent. 

While we included these corrections in our analysis we 
note that they are model dependent, and these theoretical 
uncertainties combined may account for — 5 — 10 per cent 
uncertainty over and above the statistical random errors. 


It may well be that in the future the cosmological pa¬ 
rameters will be fixed by CMB, SNe etc. Then, for fixed rea¬ 
sonable cosmological parameters, one can use redshift sur¬ 
veys to study biasing, evolution, etc. This paper is a modest 
illustration of this approach. Future work along these lines 
will include exploring non-linear biasing models (e.g. Dekel 
& Lahav 1999; Sigad, Branchini & Dekel 2001; Verde et al. 
2001) per spectral type (Madgwick et al 2002; Norberg et al. 
2002; Hawkins et al. (2001, in preparation) and the detailed 
variation of other galaxy properties with local mass density. 
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